Excitation energy after a smooth quench in a Luttinger hquid 
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Low energy physics of quasi-one-dimensional ultracold atomic gases is often described by a gapless 
Luttinger liquid (LL). It is nowadays routine to manipulate these systems by changing their param- 
eters in time but, no matter how slow the manipulation is, it must excite a gapless system. We 
study a smooth change of parameters of the LL (a smooth "quench" ) with a variable quench time 
and find that the excitation energy decays with an inverse power of the quench time. This universal 
exponent is —2 at zero temperature, and —1 for slow enough quenches at finite temperature. The 
smooth quench does not excite beyond the range of validity of the low energy LL description. 

PACS numbers: 05.30.Rt, 67.25. D-, 67.85.Hj, 64.70.Rh, 03.75.Kk 



INTRODUCTION 

Quantum physics in one dimension (ID) is much dif- 
ferent from that in higher dimensions. Many ID quan- 
tum systems, both fermionic and bosonic, have low en- 
ergy behaviour of a quantum Luttinger hquid (LL) 
No matter how complex is the underlying microscopic 
Hamiltonian, the effective quadratic LL Hamiltonian has 
only two parameters c and K, where c is a speed of its 
gapless excitations. In particular, LL description also ap- 
plies to quasi-lD ultracold al;om gases. For instance, it 
applies to a free Bose gas 0) and to bosons in a field 
of an optical lattice Q. In the case of the free Bose 
gas, the LL parameters can be related to those of the 
gas in the following way 0: K = \J c = ^^ , 

where m is the mass of a boson, p density, k com- 
pressibility, and h = 1. In particular, LL also arises 
in the system of bosons interacting via contact inter- 
actions, where [13, 13|: K = 1 + ^, c = w_f(1 — ^) . 
Here 7 = mg/h^p is interaction strength, with g being 
the contact interaction's strength, and p is linear density 
of particles. For 7 — ^ 00 we obtain a Tonks- Girardeau 
gas Q that was realized in the experiment of Kinoshita 
et al. 0. Of special interest is a Bose gas in an optical 
lattice described by a Bose-Hubbard (BH) model 



Hbh — 



Its experimental realisation has been achieved e.g. in 
Refs. When the ratio J/U in ^ is being varied, 

this model exhibits a quantum phase transition between 
the Mott insulator phase and the superfluid phase. For 
integer density we have Mott insulator pools in the p, — J 
phase diagram surrounded by the superfluid phase, see 
Figure [TJ The phase transition is of a commensurate- 
to-incommensurate type apart from situations where the 
integer density is kept fixed and the system undergoes a 
Berezinskii-Kosterlitz-Thouless (BKT) transition. These 
phase transitions can be described in terms of LL with 



K ^ 1 for commensurate - incommensurate transition 
and K = ^ for BKT transition 0, Q . The zero-gap su- 
perfluid phase can be also mapped to a LL Hamiltonian, 
and any transitions therein, driven by changing the ratio 
J /U , map to changes in the parameters c and K of the 
LL model. 
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FIG. 1: (color online). The phase diagram of a Bose- 
Hubbard model {TJ. The superfluid regime, including the 
Mott insulator-superfluid boundary, has description in terms 
of the Luttinger liquid and any trajectory in the parameter 
space of the BH model maps to a trajectory in the c — K plane 
of parameters of the Luttinger liquid. 

The ultracold gases are well isolated from their en- 
vironment and it is easy to manipulate their Hamilto- 
nian in time. Therefore, they can simulate dynamics of 
isolated quantum many-body systems driven by a time- 
dependent Hamiltonian, see Ref. [lOj for reviews. In a 
typical experiment dedicated to the dynamics, a param- 
eter in a Hamiltonian is suddenly switched between an 
initial and final value. However, even in a generic ex- 
periment, not dedicated the dynamics, a system is often 
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manipulated by continuously turning its parameter from 
an initial to a desired final value in a "quench" time tq 
(here we use the word "quench" too for historical rea- 



sons 



and assume that the state of the system is a Bogoli- 
ubov vacuum for 7fc's. In this Heisenberg picture we 
have = [a/j,iJcff]- Substituting Eq. © and using 
ll|). The idea is to prepare the system in a simple d'^^/dt = we obtain Bogoliubov-de Gennes equations 



ground state of the initial Hamiltonian and then drive 
it adiabatically to the desired ground state of the final 
Hamiltonian. Unfortunately, this ideal adiabatic quan- 
tum state preparation must fail when the system is gap- 
less either at an isolated quantum critical point or in a 
gapless phase [13,113. This failure can be quantified by 
e.g. excitation energy which is arg ued to decay with a 
power of the quench time tq |12| | , often in accordance 
with the Kibble-Zurek (KZ) mechanism An exam- 
ple of this scaling was measured recently in a smooth 
quench from the Mott insulator to superfluid 0] and in 
the Dicke quantum phase transition [8|. 

Most theoretical work has been devoted to quenches 
across a quantum critical point jTo| , where the scaling is 
determined by its critical exponents [lH . Recently, there 
has also been research on sudden interaction quenches in 
LL [13, Ell. Since most experiments begin from Bose- 
Einstein condensation, much of the (supposedly adia- 
batic) manipulation is done in the gapless superfluid 
phase of a Bose-condensed system. With this motivation 
in mind, in the following we consider a smooth quench 
of parameters in a LL which is a universal low energy 
effective model for gapless (quasi-) ID quantum systems. 

Luttinger liquid (LL). — LL Hamiltonian is 



H 



1 

r 



dx 



(2) 



where H and $ obey bosonic commutation relations. 



fc#0 



2LJ \k\ 



(3) 



*(-)--E(i) r-^-(4 + -.), (4) 

and flfe and are bosonic annihilation and creation oper- 
ators. A quench of one of the parameters in, say, the BH 
model ([U along a path in its superfluid phase, see Fig. [TJ 
maps to a quench in the LL model ([2|) along a path in the 
c — R plane of its two parameters. This map is accurate 
provided the quench does not excite high energy states 
beyond the low-energy LL model. 

In terms of and aj. the LL Hamiltonian ^ reads 



H 



= ^ c|fc 

fe#0 

R + R-^ 



R -R- 



aka^k + h.c. 



(5) 



For time-dependent R{t) and c{t) we make a Bogoliubov 
transformation 



Uk{t) Ik + w-fe(i)* life 



(6) 



dt 
C{R) 



Uk 



: c\k\ C{R) 

R + R-^ 
R - R-^ 



Uk 
Vk 

R-^ - R 

-R -R-^ 



(7) 
(8) 



Instantaneous eigenmodes of C{R) with positive norm. 



\Uk\ 



\vk\ 



{Uk,Vk) 



R + 1 R - 1 



2VR 2VR 



{U,V) (9) 



have positive instantaneous frequency c\k\. At the same 
time (V, U) is an eigenmode of C{R) with negative norm, 
|ufep — = —1, and negative frequency — c|fc|. 

Quench at zero temperature. — We drive the 
Hamiltonian ([5|) by time-dependent R{t/TQ) and 
citjTq). In the adiabatic basis ([HJ we have 



Uk 
Vk 



= ak 



-i\k\l(t) 



j\k\i(t) 



(10) 



where l{t) = dt'c{t'), and Eq. (O becomes 



ds 
ds 



a-k 



-bk e+2-/-<5l'=l^ — logi^i/2 
ds 

-ak e-2-/-«l'=l^ 4log^'/' 
ds 



(11) 



Here s = J -^r^c[t' /tq) is dimensionless time-like vari- 
able. The amplitudes ak{s),bk{s) satisfy \ak\^-\bk\^ = 1 
and initial conditions afe(— 00) = 1, bk{— 00) — 0. 

Average number of quasiparticles of momentum k ex- 
cited in the final state is Uk = |6fe(oo)p. It depends on k 
only through the product CfTQ\k\ defining a length scale 



(12) 



which is the shortest wavelength of excited phonons. 
When Uk decays with |fc| sufficiently fast, then average 
linear density of excited quasiparticles scales with tq like 



dk 
2^ 



nk 



•q 



(13) 



while the more directly measurable excitation energy 
density scales like 



'Q 



(14) 



provided that ^ A. These are universal scalings for 
quenches that do not excite beyond the range of validity 
of the LL model limited by the UV cut-off A. 



3 



ADIABATIC APPROXIMATION 



a linear ramp similar as in Ref. \14 



In the most rapid limit of a sudden quench, considered 
in Ref. 13|, we have 



lim 

CfTQ\k\- 



Uk = sinh [ log J Kf/K. 



(15) 



which is small when the relative change of K is small. 
Thus we can try an adiabatic approximation where 
|i'fcP ^ 1 and 1. Solving Eqs. (fTTj) perturbatively 

to leading order in bk yields a Fourier transform 



nk 



ds 



(16) 



This rik is small and the adiabatic approximation is self- 
consistent when relative changes of K during a quench 
are small. Moreover, even when they are large Eq. (fT6|) 
is still accurate for large enough ^|fc| where Uk is small, 
see the following examples and Figure [2] 



tanh 
A = 0.1 
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FIG. 2: (color online). The figure displays momentum de- 
pendence of excitation density Uk as a function of 2^|fc| for 
the tanh quench (left panels) in Eq. (|17p and the linear ramp 
(right panels) in Eq. (|18|l . The upper row shows results for 
small quenches with ^4 = 0.1 when the adiabatic approxima- 
tion (|16|l agrees with numerical evaluation of equations (lll|l . 
For a greater A = 0.8 (lower row) it remains valid for 2^|fc| 
large enough, but not for small 2^|fc| where it underestimates 
Uk that tends to (|15|l for 2^|fc| — >■ instead of log^ ^jKflKi 
from a perturbative Eq. (I16p . The insets show corresponding 
plots in a linear scale. 



Here we consider three examples of quenches around a 
K^: a smooth tanh quench 



\o%K^^\s)IKq = A tanh(s) , 



(17) 



logX(2)(5)/^g 



-1 , when s < —1, 
A { s , when — 1 < s < 1, 
1 , when s > 1, 



and a smooth shake beginning and ending at Kq 

A 



cosh(s) 



(18) 
(19) 



Corresponding power spectra (IT6t are 



,(1) - 



,(2) 



,(3) 



A^ 
A^ 
A^ 



sinh (7r^|fc|) 
sin {2^\k\ 

<\k\ 
cosh (7r^|fc|) 



(20) 



They are exponentially localized on a scale fc cx ^ ^ ex- 
cept the linear ramp (2) whose discontinuous time deriva- 



tive results in a fat high energy tail n 
linear density of excitations (|13l) reads 



(2) 



A2 
6^' 



4^ 



,(3) 



12^' 



The 



(21) 



when ^ 1. All three densities scale like Tq^ in agree- 
ment with the nonperturbative Eq. (1131) . 
The excitation energy density (IT4)) reads 



.(1) 



r(2) _ 



3C(3) A^ 
27r3 ^2 
J_ A2 
8^ ^ 



.(3) 



2 log(2^A) 



9C(3) A2 
87r3 ^2 



(22) 
(23) 



when 5 A 3> 1. Here C, is the Riemann zeta function. The 
energy scales like Tq^, in agreement with the nonpertur- 
bative Eq. (fT4| . except case (2) when it logarithmically 
diverges with the cut-off. The discontinuous time deriva- 
tive of the linear ramp (1181) is probing non- universal high 
energy excitations beyond the range of validity of the ef- 
fective low-energy LL ([2]). For the linear ramp the LL 
is not a self-consistent approximation and the scaling of 
the excitation energy with the ramp time tq, if any, can 
be different from the universal Tq^. A non- universal ex- 
ponent was observed in numerical simulations of linear 
ramps in the BH model 14| . 



CORRELATIONS 

The excitation does not change the quasi-long-range 
character of correlation functions. For example, after a 
quench is completed a correlation function for density 
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fluctuations smeared over distances longer than interpar- 
ticle distance 5p{x) = -dx^{x) reads 



C{R) 



Kt 



{5p{x 

\k\dk 
47r3 



+ R)5p{x)) 
[1 + 2nfc 



2v^cos(^fc(i)] ,(24) 



where Lpk{t) = 2c/i|fc| + (/3"(^|fc|). In the ground state, 
when rife = 0, its tail decays like Cf{R) 



1/A. The excitations Uk > add a correction 



for 



G 



2c/t - R 



G 



2cft + R 



(25) 

Here F{z) and G{z) are real functions originating from 
the Uk and y/nk terms in Eq. ([24]) respectively . 

There is analogy to the quasiparticle horizon effect ■ 
The G-terms in ((25|) describe shock waves that originate 
from correlated pairs of quasiparticles with momenta ±fc 
excited during a quench whose separation grows like 2c ft. 
Their width ~ ^ is the shortest length on which quasipar- 
ticles excited in time tq can be localized in space. The 
_F-term in is an additive correction to Cf{R) that 
remains inside the quasiparticle horizon, R ^ 2c ft, after 
the shock waves go away. When no > in (|15p. then a 
tail of the remaining correlator is C{R) — (l + 2no)G/(i?) 
for R ^ S,. The R^^ tail of the ground state correlator 
is amplified by a factor (1 + 2no) — ^{Kf/Ki + Ki/Kf). 
When Tio = in ([TS]), as for e.g. the shake (|19I20I) . then 
the F term decays faster than and the tail of C{R) 
is the same as in the ground state: C{R) = Cf{R) for 
C<i?<2c/t. 

Despite all its interesting physics, the additive correc- 
tion ((25|) does not cut the algebraic i?~^-tail of the corre- 
lator in the quasi-long-range ordered ground state. This 
contrasts with the KZ mechanism where correlations af- 
ter a quench from a disordered to an ordered phase decay 
exponentially [Ift] . 



QUENCH AT FINITE TEMPERATURE 

Since real experiments are done at finite temperature, 
we generalize to finite T, where the excitation spectrum 
reads 



nk{T) = 



n, (1 + 2nr) 



(26) 



Here rik — |6fc(oo)p is the excitation spectrum at T = 
0, and nf^ = [exp{ci\k\/T) — 1]^^ is the initial thermal 
distribution which can be also reinterpreted as a final 
thermal distribution n^^ = [exp(c/|fc|/T/) — 1] ^ with 



Tf = -^T 
c, 



(27) 



This temperature rescaling is the only effect in the adia- 
batic limit tq — oo when — ^ 0. 



Here we are more interested in the non-adiabatic exci- 
tation above this thermal background: 



nk{T) 



RE 



nk 



(l + 2nn 



(28) 



When compared with the bare Uk at T = 0, it is ampli- 
fied by a Bose enhancement factor 1 + 2n^^ making the 
excitation more significant than at T = 0. For relatively 
fast quenches, with tqT <^ 1, the distribution nk extends 
up to |fc| ~ where the enhancement 1 + 2n^^ « 1. 
Consequently, the non-adiabatic excitation energy den- 
sity e{T) is roughly same as at T = 0. In contrast, for 
relatively slow quenches with 



TQ T » 1 



(29) 



the fc — > singularity of the Bose enhancement factor 
has more qualitative consequences. In this regime the 
non-adiabatic excitation energy density is 

eiT) = f^^Cf\k\nk{l + 2nr) 



/^§C,|fc| Uk 



2T 



f 



Cf\k\ 



= 2Tf ne, , (30) 



where riox is the density of excited quasiparticles at T = 
in (|13p . The energy is not only larger than its T = 
counterpart (|14p . but its decay with tq is also less steep: 

e{T) ~ r-i (31) 
from (|13I30P instead of the r'^ at T = in (Hi)). 

CONCLUSION 

We derived universal dynamical exponents for the scal- 
ing of excitation energy with a quench time when the 
quench is slow and smooth enough not to excite beyond 
the quadratic Luttinger liquid model. Due to the singu- 
larity of the Bose enhancement factor, the exponents are 
different at zero and finite temperature. In both regimes 
they can be used to test the low energy LL description 
of quasi- ID ultracold atomic systems. 
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